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ABSTRACT 


V 

This  paper  solves  a  targeting  problem  exactly  end 
compares  the  exact  solution  with  several  reasonable  expected- 
value  approaches  in  order  to  gauge  their  accuracy.  The 
scenario  consists  of  air-to-surface  missiles  fired  simul¬ 
taneously  at  a  ship  formation  defended  with  point  defenses. 
The  problem  is  to  find  the  probability  that  a  given  ship 
will  survive  the  raid.  Each  ship  is  present  with  a  given 
probability  and  possesses  a  targeting  "weight"  determining 
the  likelihood  that  it  will  be  selected  for  targeting  among 
the  other  ships  present.  The  method  can  be  extended  to 
further  layers  of  defense  pnder  certain  conditions.  Each 
solution  is  presented  in  2  forms:  a  finite  sum  and  a  definite 
integral,  each  being  easier  to  handle  than  the  other  under 
certain  conditions.  Examples  are  given.  One  important  con¬ 
clusion  is  the  spectacular  breakdown  of  reasonable  expected- 
value  approaches. 
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A  TARGETING  PROBLEM:  EXACT  VS.  EXPECTED-VALUE  APPROACHES 

INTRODUCTION 

The  purpose  of  this  paper  is  twofold:  first,  to  solve  a  targeting 
problem  exactly,  and  second,  to  compare  the  exact  solution  with  several 
possible  expected-value  approximations  in  order  to  gauge  their  accuracy. 

The  targeting  problem  consists  of  a  wave  of  bomber  aircraft  simul¬ 
taneously  firing  missiles  at  a  ship  formation.  )}efore  the  first  bomber 

'  i 

raid,  all  the  ships  are  present,  bqt  after  the  first  raid  a  ship  is  present 
only  with  a  certain  probability  — -  depending  on  whether  it  was  sunk  or  not. 
Strictly  speaking,  these  probabilities  are  not  independent  (a  missile  can  hit 
only  1  ship:  when  it  does,  it  is  no  longer  available  to  hit  another  ship). 
However,  if  they  are  assumed  independent,  multiple  raids  can  be  analyzed 
sequentially  with  the  method  presented  here.  This  method  exactly  calculates 
the  probability  that  a  given  ship  is  present  after  a  raid  from  the  probabili¬ 
ties  of  ship  presence  before  the  raid  and  other  relevant  parameters  (number 
of  missiles  launched,  ship  targeting  weights,  and  ship  vulnerability  data). 

The  bombers  typically  face  many  layers  of  defenses,  but  the  model 
described  here  considers  only  those  missiles  facing  the  last  layer,  namely 
the  ships’  point  defenses,  if  any.  Nevertheless,  the  method  can  be  extended 
to  further  layers  ship  defense  --  such  as  SAM  barriers  or  outer  air  battle 
of  fighters  vs.  bombers. 

< 

Six  expected -value  approaches  are  presented  as  possible  alternatives 
to  the  "exact”  method.  Their  obvious  advantage  is  that  they  are  easy  to  use 
and  require  far  less  computer  time  and  space.  Two  numerical  examples  help 
determine  which,  if  any,  of  the  6  approximations;  are  acceptable  and  under  wha£ 
conditions  they  are  acceptable. 


THE  PROBLEM 

The  object  is: 

(1)  To  calculate  the  probability  chat  a  given  ship  in  a  formation 
is  targeted  by  a  given  number  of  incoming  antiship  missiles  (ASMs) ,  and 

(2)  To  calculate  the  probability  that  a  given  ship  will  survive 
the  raid. 

(We  shall  see  that  it  is  not  always  necessary  to  calculate  (1)  in  order 

to  obtain  (2).)  The  inputs  are: 

n  :  maximum  number  of  ships  in  formation 

a.  :  probability  that  ith  ship  is  present  in  formation 
(the  a^s  are  assumed  independent) 

w^  :  weight  of  1th  ship  (a  positive  integer) 

M  :  total  number  of  ASMs  launched  (the  missiles  are 
•  considered  independent) 

L. (m) :  lethality  function,  i.e.,  probability  that  the  ith  ship 
is  sunk  if  hit  by  ra  missiles  (or,  if  m  missiles  attempt  to 
penetrate  its  point  defenses). 

The  constant  weight  assigned  to  each  ship  determines  the  probability  that  the 
ship  will  be  targeted.  For  example,  the  probability  that  the  5th  ship 
is  targeted  given  that  the  2nd,  5th,  7th,  and  10th  ships  alone  'are  present  is 
v5/ (w2  +  w5  +  w7  +  w10)  . 
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THE  SOLUTION 


All  solutions  will  be  expressed  both  as  suns  and  integrals.  In 
general,  the  larger  the  number  of  ships  n  the  longer  it  takes  to  calculate 
the  sum  (2°^  terras)  and  the  larger  the  number  of  missiles  M  the  longer  it  takes 
to  calculate  the  integral  (up  to  M  terns  in  the  integrand) . 


Targeting  Probabilities 

P  (ra|M)  =  probability  that  ith  ship  is  targeted  t>y  exactly  m  missiles 
given  that  M  missiles  are  launched  and  given  that  ship  i 
is  present 


3  y  *  ^  1  a ,  .a .  •  *  •  a , 

k=l  *  {!,..» ,n;excppt  i}  J2  J3  Jk 

•  (1-a.  ) . . .  (l~a.  )(”)p“  (l-P.)”"0  .  (1) 

jk+l  Jn  *  1  -  1 


where 


P 


i 


w.+w. 
1  J 


(2) 


This  sum  can  be  expressed  ps  the  following  definite  integral  for 
m/0:  ^ 

Pi(m|M)  -  wi(”)  /  fi(x)dx  •  (3) 

0 
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where 


w._i  “  w. 

=  g^(x)x  1  n  (a  x  j  +  1-a  ) 

1  1  j-l;*i  j  j 


,M-mw  .,4  (-wiL°gx) 
8i(x)  -  L  <  i  )(_1)  " '(Jl+m-l)  * 

Jl=0 


4+m-l 


For  m^O, 
P 


.(0|M)  =  1  +  w^  f dx  h.(x)x  1  II  (a  x  J  +  1-a.) 


(4) 


(5) 


(6) 


where 

M 

h,(x)  =  52  (^)(-l)A(-wiLogx)A“1/(fc-l)!  (6a) 

It  is  observed  that  the  distribution  is  unimodal  (rises,  peaks,  then 
falls),  except  for  the  fact  that  it  may  rise  abruptly  at  the  endpoint  m=M. 
This  corresponds  to  the  configuration  where  ship  i  is  the  only  one  present, 
in  which  case  it  catches  all  the  missiles. 


Survival  Probabilities 

The  probability  that  ship  i  will  be  present  after  the  raid  is: 


a^  -  af[l  -  Pk1(M)]  ,  (7) 

where 

PKi(M)  .=  probability  that  ship  i  is  sunk  given  that  it  is  present 
and  that  M  missiles  are  launched  on  the  formation 

M 

-  £  L^m)  P^mlM)  .  (8) 

m**0 
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Note  that,  strictly  speaking,  this  scheme  cannot  be  iterated  for  several 

missile  raids  because  the  a'  are  pot  independent  (as  the  a  were 

1  1  * 
assumed  to  be)  . 

Specific  forms  for  the  lethality  function  L^ra)  often  enable 
one  to  perform  the  above  sum  directly,  without  having  to  specifically 
calculate  the  targeting  distribution  Pi(m|M).  Consider  the  following  form: 


L^m)  -  1  -  (l-d1)m  +  cim(l-b1)m  ,  (9) 

where  <1^  and  b^  are  free  parameters  determined  by  curve-fitting  the 
lethality  data  (frequently,  a  simple  power  law,  i.e.,  c^O.is  sufficient. 
Then 


Pv .  (M)  *»  22  a4  (l-a.  )••• 

k=l  “  (l, . . . ,n;except  i}  J2  Jk  Jk+1 


(1"aj  ^  f1.  “  <1-diPi>M  +  Mc1p1U-bi)(l-b1Pi)M"1j 


with  pi  given  in  equation  (2).  The  integral  form  is 


(10) 


FK1CM) 


-  v  f  dx  u.Wx"1'1  8  ‘v"3  + 

0 J  ‘ 


(11) 


To  see  this,  consider  a  formation  of  2  ships  with  equal  weights  and 
ai=(l,l)  targeted  by  1  missile  with  lethality  1  (nuclear  missile). 

Then  .5)  •  If  the  a^  were  independent,  the  probability  that 

both  ships  would  be  sunk  f ( 1— . 5) ^  *  .25]  would  be  nonzero,  clearly  an 
impossibility  since  the  missile  can  only  hit  one  ship. 

The  results  shown  here  are  strictly  correct  only  if  the  initial  probabilities 
of  presence  (a^)  are  truly  independent.  This  means  that  the  a^s  cannot 
themselves  be  the  pesuit  of  a  missile  raid.  They  could  be,  for  example,  the 
result  of  independent  attacks  on  each  ship  prior  to  the  bombing  raid. 
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where 


If  the  lethality  data  are  too  irregular  to  be  fitted  by  (9),  other 
formulas,  with  more  free  parameters,  can  be  constructed  which  might  enable 
one  to  carry  out  the  sum  in  (8)  exactly.  These  lethality  formulas  must 
satisfy  <XL^(m)<l,  1^(0)  =  0,  and  L.(»)  =»  1.  An  example  is 

■L±(m)  -  1  -  (l-d1)in  +  P1(m)(l-bi)m  +  P2(rn)(l-b’i)m  +  ..., 


coefficients  are  free  parameters  depending  on  i. 


Even  when  specific  forms  for  the  lethality  function  L^fm) 
allow  us  to  bypass  the  calculation  of  the  targeting  distribution  P^mlM), 
some  useful  statistical  information  can  still  be  obtained.  For  example: 


*  m^  =  expected  number  of  missiles  targeting  ship  i  given 

that  it  is  present 

=  MPi(l|l)  .  (14) 

Calculation  of  P.(ljl)=p.  is  particularly  simple 
(reference  1). 

•  The  expected  number  of  missile  hits  on  6hip  i  can  also 
be  calculated  without  knowing  _P.(mjM).  If  there  are  no 
point  defenses,  it  is  simply  m.1.  If  there  are  poirjt 
defenses  (which  wquld  be  included  in  the  lethality  functions) , 
it  is  obtained  by  matching  the  probability  of  sink  P^ 

with  the  curve  giving  probability  of  sink  versus  number  q f 
hits  received. 


PROOFS 

Proof  of  Targeting  Distribution  Formulas 

To  prove  the  6um  result  in  equation  (1),  consider  a  subset  of  k  ships 

containing  the  ith  ship.  The  ships  in  the  subset  are  labelled  i» » • *  * 

The  probability  that  this  subset  is  realized  is  a. a.  ...a,  (1-a,  )...(l-a  ), 

1  J2  3k  3k+l  Jn 

and  the  probability  that  a  given  missile  will  target  the  ith  ship  from  this 

subset  is  p,  3  w./(w,+w,  +. . .+w  ).  Since  the  missiles  are  independent,  the 

1  1  1  ^2  ^k 

distribution  for  m  missiles  impacting  is  binomial  with  probability  p^  .  One 
must  then  sum  over  all  possible  ways  of  forming  a  subset  of  k  ships  containing 
the  1th  ship,  and  sum  also  over  all  possible  sizes  of  the  subset  (k  varying  from 
1  to  n) .  Each  of  the  2n  *  possible  configurations  of  ships  —  from  the 
ith  ship  alone  to  all  ships  present  —  will  then  have  been  included. 
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Let  us  now  turn  to  the  integral  representation  (3).  The  case  for 
1  missile  (m=M=l)  was  treated  in  reference  1,  and  the  procedure  is  generalized 
here.  We  can  establish  the  following  lemma. 


Lemma 


The  following  sum  : 


^  S  a-t  •  ■  .a .  (l~a .  )  . . .  (1-a .  )$(P. ) 

k=l  { j ~ , . . .  ,j.  }={  1, ...  ,n;except  i}  2  k  k+1  3n 

1  (15) 


where  \p  Is  an  arbitrary  function  and  p.  =  w./(w.+w.  +...-K?.  ),  is  equal 

1  1  1  3  2  3k 

to  the  following  integral: 


w.  f  dx  F  (x)x  1  II  (a  x  3+l-a  ),  (16) 

1  J  1  j«l;*i  3  3 

0 


where  F^(x)  is  the  inverse  Laplace  transform  of  i^(l/s)  evaluated  at 

t  ®  -w^Log  when  such  a  transform  exists.3--  [The  Laplace' transform  of 
00 

q(t)  is  J  e-Stq(t)dt]. 

0 


F^(x)  can  be  a  distribution,  rather  than  a  function  defined  at  every  point. 
All  that  is  required  is  that  the  Integral  in  (16)  be  defined.  See  below 
for  an  example  where  F  is  the  distribution  called  the  delta  "function". 
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terms 


Proof  of  Lemma.  Both  the  Integral  and  the  sui,.  have  ^ 


A  typical  term  in  the  integrand  of  (16)  is 

w.-fw.  +...+W,  -1 

1  d  2 

w.a.  ...a.  (1-a,  )...(l-a,  )  x  F  (x) 

32  3k  3k+l  3n  1 


To  match  this  with  the  typical  term  in  the  sum  (15),  we  must  have 


“i  / 


1  W.+W.  +...+W,  -1 

-  i  1 


Jk  Ft(x)d ,  , 

\  1  h  K' 


(17) 


i 

If  we  make  the  change  oj;  variable  x  *  y,  then  define 

1/w. 

s  =  1  +  (w.  +...+W  )/w  And  <P.  (y)  ■  F .  (y  x),  the  problem  reduces  to 

32  3k  1  1  1 

finding  a  function  4>^  independent  of  s  such  that 


1 

/  ys-1  <J>1(y)dy  -  Mj) 
o 


With  the  change  of  variable  y*e"t,  this  becomes: 

/» 

e'st<f>(e't)dt-4)(i)  . 

i.  ® 

0 

The  left-hand  side  is  the  definition  of  the  Laplace  transform  of  (^(e  C), 

wi 

Substituting  beck  t  “  -Log  y  and  y  *  x  gives  the  result  stated  in  the 
lemma. 


In  this  case,  we  have 
4><P i>  *  P  ”  (l-Pi)”-”  • 
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so 


*$>  ■  ^ 
s 


M-m 


M-m 

E,M-m.  ,  , .  I ,  m+£ 
(  p  )(-l)  Is 

£= 0 


(18) 


Taking  the  inverse  Laplace  transform  of  the  right-hand  side  gives 
M-m 

Me“C>  -  Z  (MIm)(-l)£  t<l'hn‘1/(£-Hn-l)!  (19) 

£=0 

wi 

for  m^O.  Substituting  back  t  =  -Log  y  and  y  =  x  gives  the  g^(x) 

in  equation  (5).  When  m=0,  the  £=0  term  in  (19)  is  undefined. 

Equation  (18)  shows  that  this  term  is  the  inverse  Laplace  transform  of  1, 

which  the  delta  function^  6(t)  .  Thus, 

M 

^(e"*)  -  6 (t)  +  £  (^)(-l)£  t*-1/(£-l)!  (20) 

£=1 

for  m=0.  Substituting  back  to  y,  then  x,  and  putting  the  resulting  g^x) 

in  (4)  yields  the  expression  (6)  for  P  (0|M)  .  A  change  of  variable 

b  1 

-Log  x  =  u  and  the  fact  that  J  5(x)f (x)dx  -  f(0)  if  0efa,b]  must 
be  used. 


Note  that  gi  is  essentially  a  Laguerre  polynomial.  Indeed,  the 
Laplace  transform  of  tM-1(M-m)!  L^“*  (t)/(K-l)!  is  (s-l)M-m/sM,  which 
is  the  right-hand  side  of  (18)  (reference  2,  p.  191").  Thus,  for  m^O  > 

g±M  -  (-w.Log  x)M_1  (M-m)!  Ljj^  (~wi  Log  x)/(M-l)l  (21) 


The  step  function  9(t-k),  which  is  1  for  t>k  and  0  otherwise,  has  the 
Laplace  transform  e”^s/s.  Its  derivative  with  respect  to  k,  -5(t-k),  has 
the  Laplace  transform  -e“^s.  As  k-*-0+,  6(t)  has  the  Laplace  transform  1. 
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where  L  is  the  Laguerre  polynomial. 


l*<*>  $  E  O (-*>■/*: 

n  ~  n-m 

m=0 


Proof  of  Survival  Probabilities  Formula 

- - - - - - - - - - - f - - - 

For  a  lethality  function  L^(ra)  given  by  (9),  one.  can  readily 
establish  that  P^{M)  is  given  by  (10)  by  using  the  following  results: 


V'  m  .M-m  , 

m=0 

M  % 

E.M..  mm..  .M-m  .  .  ,  .M 

<m)  ^  P±  (1-Pj)  -  (c^  +  1  -  Pt) 

m=0  1 

M 

S.M.  m  m  .M-m  ,  .  M—  1 

(m)ma  pt  (1-P1)  *  Motpi(api  +  1  -  Pi> 


(22) 


(23) 


(24) 


Equations  (22)  and  (23)  are  simple  binomial  expansions  of  (a+b)n. 
Equation  f24)  can  be  obtained  from  (23)  by  differentiating  with  respect 
to  a  then  multiplying  by  a. 


To  find  the  u^(x)  in  the  integral  form,  we  proceed  as  before, 
reference  to  the  lemma  established  earlier,  we  have: 

*(Pt)  "  1  "  (1-diPi)”  +  MciPi(l-b1)(l-biPi)M“1 


With 


(25) 


so 


/s-d  \M  (s-b.) 

^  [ _ 7;  I  v».  /i  l  \  1 

Vs/ 


M-l 


+  MCi(l-bt) 


M 

»  M  (-d,)*  M-l  „  ,  (-b. )*■ 

-f  r~  +  E  <Y)  -TTi- 

£-1  s  •  £»0  s 


(26) 


Taking  inverse  Laplace  transforms  as  before  gives  <p^(e  t),  and  hence 
u^(x)  .  A  change  of  index  in  the  first  sum  then  puts  it  in  the  form 
shown  in  (12). 


Formula  (14)  is  proved  by  observing  that 


X)  a-Pirm  =  MPi  (i.e.,  (24)  with  a-1) 

m=0 


PROPERTIES  OF  THE  INTEGRAL  REPRESENTATIONS 

The  integrand  f.(x)  in  (3)  has  the  following  properties: 

1  Wj-1 

(1)  fi(0)  *  0  if  w^>l  because  of  the  presence  of  x  and  the 

N 

fact  that  for  any  positive  integer  N,[(Log  x)  x] tends  to  0 
when  x  tends  to  0.  If  w^l,  then  f^O)  can  be  infinite. 
Since  this  should  be  avoided  in  numerical  integration,  it  is 
preferable  to  take  all  weights  w^  larger  than  1,  which  does 
not  change  the  end  result. 


<2>  Et(1>  -{?  “  t\ 


(3)  The  polynomial  in  x  entering  the  expression  of  f^(x)  has 

positive  coefficients,  and  hence  no  roots  for  x  in  [0,1].  The 
function  g^(x),  which  is  a  polynomial  in  the  positive  quantity 
y=  -  Log  x  whose  coefficients  have  alternating  signs,  has 

M-m  positive  roots  y  (reference  2,  p.  204). 

»  t  __  —  —  - 


Therefore, 


f^(x)  has  M-m  roots  (xf  *  e  )  between 


0  and  1,  excluding  0  and  1.  Thus,  the  graph  of  f^(x)  ia 
oscillatory  with  (M-tu)/2  cycles.  The  graphs  of  hj,(x)  in 
(6)  and  u^(x)  in  (12)  are  also  oscillatory. 
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(4)  The  case  M=m=l  is  particularly  simple  to  integrate 

numerically  and  was  treated  in  reference  1.  f(x)  is  a  polynomial 
with  positive  coefficients,  and  hence  its  graph  is  concave 
upwards.  The  maximum  error  using  the  trapezoidal  rule  or 
better  was  given  as  [ 2f (0)  +  f(l)  -  f (1-1/N)]/2N,  where  N  is 
the  number  of  (equal)  intervals  used.  A  refinement^",  due 
to  David  A,  Perin,  gives  a  maximum  error  of  only 

E  =  [f CO)  +  f(l)  -  f (1-1/4N)  -  f (1/4N) ] / 2N  .  (27) 


Verifications 

Any  program  which  calculates  P^(mjM)  should  include  the  following 
simple  checks: 


M 

•  £  Pj.  (rajM)  -  i 

m=0 

M 

•  £  m  M*|M)  -  M  P.  (1 1 1) 

m=0  1  1 


(28) 

(29) 


The  quantity 


P1(l|l)  =  Pj  must  be  calculated  separately  as  described  earlier. 


Instead  of  extending  the  chord  from  x^  to  t0  obtain  a  lower  bound 

on  the  next  interval,  as  was  done  in  reference  1,  consider  the  tangent  to 
the  curve  at  each  interval  point  xj  (except  the  endpoints  0  and  1),  and 
extend  it  until  it  crosses  the  raid-interval  lines  on  either  side 
[y  “  (x±  +  x±„i)/2  and  y  ■  (x^  +  Xi+^)/2]  .  The  area  of  the  resulting 
right  trapezoid  is  siraply  f(xi)/N.  This  gives  a  lover  bound  for  all  but 
2  half-intervals,  one  at  each  endpoint.  Drawing  the  tangent  at  the  mid¬ 
point  of  each  half-interval  [at  Xj_/4  and  (3  +  X}j_i)/4]  gives  the  total 
lower  bound: 

/ fWdx  >  £[f(£)  +  ...  +  fCl-^)]  +  ^  tf(^)  +  . 

0 
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EXPECTED-VALUE  APPROACHES 

Expected -value  constructions  of  the  targeting  distribution 
P^mjM)  consist  in  calculating  an  average  probability  p^  that  a  single 
missile  will  target  ship  i  given  that  ship  i  is  present  and  then  constructing 
a  binomial  distribution  based  on  this  average  probability.  Thus 


P1(»|M)  =  <”>  p”  d-Pi) 


,M-m 


(30) 


Some  possibilities  for  p^  are: 


Expected-Value  Approach  No.  1:  (Semi-Exact) 


> -P  (l|l)-  £ 

A  ^  * _ 1 


2-# 

k=l  =  {1* . . . »n;except  i}  J2  Jk 


w, 


(1-a.  ) . . , (1-a .  ) 


'k+1 


j  w.+w.  +...+W, 

n  1  U  1 


w.,-1  n  w. 

(ax  J+l-a,)dx 


1  u  n 

w  f  x  1  n 
1 J  j-i,i 


(31) 


Expected-Value  Approach  No.  2: 


aiwi 


i  ^n 

£  a1wi 

J"1  3  3 


(32) 
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No.  1  is  referred  to  as  "serai-exact"  because  the  it  uses  is 

exact.  No.  2  takes  the  ratio  of  the  expected  weight  of  ship  i  to  the  total 

expected  weight  of  the  formation.  No.  3  takes  the  full  weight  of  the 

ith  ship  (since  p^  assumes  ship  i  is  present)  divided  by  the  expected 

weight,  and  normalizes  the  result,  so  that  the  probability  of  targeting 
n 

something,  namely  £  ai  *  *s  e<lual  t0  the  probability  that  something 

i=l  _ 

is  present,  namely  1  -  II(l-a  ).  Its  drawback  is  that  p  can  be  larger 

s  s  1  ' 

than  1  (if  the  a^s  are  small  enough).  No.  4  argyes  that  since  the  ith 
ship  is- assumed  present  in  p^,  it  should  contribute  Its  full  weight  when 
the  expected  weight  of  the  formation  is  computed  in  the  denominator.  No.  5 
normalizes.  No.  4.  No.  6  assumes  that  all  ships  are  present  for  purposes 
of  calculating  the  conditional  probability,  with  a  view  of  correcting  this 
assumption  when  all  P^OnlM)  are  multiplied  by  a ^  at  the  end  to  yield 
the  unconditional  probability.  All  approaches,  exact  or  expected-value, 
must  reduce  to  No.  6  when  all  a^s  are  indeed  1.  No.  6  gauges  the  sensitivity 
of  the  resulting  probability  of  targeting  ship  i  given  that  it  is  present 
to  whether  or  not  ship  i  is,  in  fact,  present. 

Critique  of  the  Expected-Value  Approaches;  Examples 

We  will  compare  the  various  approaches  in  2  examples.  Both 
include  n=ll  ships.  In  the  first  the  formation  is  attacked 
with  M =  12  missiles.  In  the  second  the  weights  are  slightly 
changed  and  M  =  40  missiles  attack.  The  data  are: 
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»a  I o. 


a±  -  .379,  .249,  .605,  .553,  .458,  .6?5,  .258,  .695,  .437,  .460,  .553 


Vi3  2,  2,  1,  1,  1  ,  1,  1,  1,  1,  1,  1 

M  =  12 


Example  2: 

a±  =  .379,  .249,  .605,  .553,  .458,  .695,  .258,  .695,  .437,  .460,  .553 
wi-  4,  4,  1,  1,  1,  1,  1,  1,  1,  1,  1 
M  -  40 


The  exact  and  tbe  6  expected-value  targeting  distributions  for  ship  1, 
P^(m|M),  are  shown  in  figures  1  and  2*.  We  can  make  the  following  observations: 

•  When  only  12  missiles  are  launched,  many  expected- 
value  approaches  come  reasonably  close  to  the 
correct  result.  When  as  many  as  40  missiles  are 
launched,  all  expected-value  approaches  break  down. 

They  all  lead  to  narrow,  sharply-peaked  distributions 
whereas  the  correct  distribution  is  noticeably  spread  out. 


^The  p ;  on  which  the  binomial  distributions  are  based  are: 
Expected^value 

method  number  '  1 _ 2 _ 3 _ 4  5 _ 6 

L  (example  1)  .295  .127  .335  .277  .306  .154 

.  (example  2)  .435  .210  .553  .412  .474  .235 
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FIG.  1:  Targeting  Distribution  In  Example 


Targeting  Distribution  in  Example 


Approaches  Nos-  2-  and  6  lead  very  far  from  the  correct 
result  in  both  cases. 


•  The  serai-exact  result  (No.  1)  is  closest  to  the 
exact  result  in  both  cases. 

•  No.  3  peaks  later  than  all  the  others. 

•  Nos.  1,  4,  and  5  are  close  to  one  another  in  both  cases. 


The  wider  spread  of  the  exact  distribution  can  be  predicted  by  cal¬ 
culating  variances.  For  example,  the  variance  of  the  exact  distribution 
(calculated  by  James  K.  Tyson)  is: 


where 


a2  =  m  p. (1-p.)  +  M(H-l)a  2 
in  _ l  l  p 

a  2  5  p,2  -  p  2 


P  - 1  ‘i 

exact  distribution,  which  is  a  single  binomial,  is 


is  independent  of  M,  but  the  variance  of  the  semi- 


m 


M  P1(l“Pi) 


Thus,  the  exact  distribution  is  always  more  spread  out  than  the  serai-exact 
expected-value  distribution,  and  the  difference  between  the  2  standard 
deviations  grows  as  M,  the  total  number  of  missiles.  For  a  very  large 


number  of  missiles,  one  would  witness  a  spectacular  breakdown  of  all  the 
expected-value  approaches  considered  here. 

Let  us  now  introduce  the  lethality  (Equation  (8))  to  find  the 

probability  P  ,  that  ship  1  is  sunk  given  it  is  present.  Two 
K-L  ^ 

lethality  functions  are  considered:  L^(m)  =  1-(1-  .01)  and 
L^(m)  =  l-(l-.03)m.  The  results  for  example  2  are  given  in  the 
table  below. 


PK1 

Lethality  function 

I*,  (m) 

Exact 

No.  1 

(semi- 

exact) 

No.  2  No.  3  No.  4  No.  5 

No.  6 

- -  l 

1  -  <l-.01)m 

.159 

.160 

.081  .199  .152  .173 

.090 

1  -  (l-.03)m 

.404 

.408 

.223  .488  .392  .436 

.247 

As  can  be  seen,  the  serai-exact  result  is  extremely  close  to  the  exact 

result,  and  Nos.  4  and  5  are  not  far  behind.  The  other  expected-value 

methods  break  down.  It  appears  therefore,  that,  in  this  particular  case, 

expected-value  distributions  having  the  same  mean  as  the  exact  one  yield 

results  close  to  the  exact  one  for  power-law  lethality  functions,  however, 

this  is  not  true  for  other  realistic  lethality  functions.  For  example, 

consider  the  case  of  the  cookie-cutter  lethality  function, 

L^(ra)  =  0  for  an^  1  for  m>mQ,  corresponding  to  a  saturation  effect 

with  nuclear  missiles  (the  point  defenses  destroy  any  number  of  missiles 

smaller  than  mo,  but  if  more  than  raQ  engage,  the  excess  will  hit  the 

ship  and  sink  it  with  cartainpy) .  The  sink  probability  is  then  simply 
M 

P  (M)  =  ^2  P » (m | M) ,  which  can  vary  widely  depending  on  m  and  the 

m=nio  X  ° 

approach  used,  as  figures  1  and  2  testify. 

The  reason  why  the  sink  probabilities  given  by  some  expected-value 
results  are  close  the  exact  sink  probability  in  this  case  is  that  in  the 
interval  where  most  of  the  distribution  is  concentrated  (m  between  12  to  24) 
the  lethality  function  does  not  vary  dramatically,  so  that  the  shape  of 
the  distribution  is  not  a  dominant  consideration. 


Conclusion  on  Expected-Value  Approaches 

The  conclusions  to  be  drawn  from  the  above  are: 

•  When  the  object  is  to  calculate  the  targeting  distribution, 
all  expected-value  approaches  fail  significantly  for  large 
numbers  of  missiles  ("40),  but  some  are  successful  for 
small  numbers  of  missiles  (~12). 
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•  When  Che  object  is  to  calculate  the  probability  of  sinking 
a  ship,  all  expected-value  approaches  fail  significantly 
for  some  lethality  functions  (such  as  a  cookie-cutter), 
but  some  are  successful  for  other  lethality  functions 
(such  as  power  laws). 

#  The  semi-exact  expected-value  approach  (No.  1)  yields 
results  closest  to  the  exact  approach.  Approaches  Nos.  4 
and  5  are  not  far  behind.  The  other  approaches  fail 
significantly  in  all  cases. 

As  can  be  seen,  results  yielded  by  reasonable  expect'ed-value  approaches 
can  show  large  errors  when  compared  to  results  yielded  by  the  corresponding 
exact  approach.  This  is  all  the  more  sobering  as  in  many  cases  the  problem 
is  so  complex  (e.g.,  a  full  campaign)  that  an  exact  approach  is  almost 
impossible  to  formulate,  let  alone  to  use  for  calculations. 

More  details  and  computer  programs  can  be  found  in  reference  3. 
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